4th. November 2009 11:12 Ratcliffe 



NON-UNIVERSALITY AND EVOLUTION 
OF THE SIVERS FUNCTION 



PHILIP G. RATCLIFFE 

Dipartimento di Fisica e Matematica, Universitd degli Studi dell'Insubria, Como 

and 

Istituto Nazionale di Fisica Nucleare, Sezione di Milano-Bicocca 
OLEG TERYAEV 

Bogoliubov Laboratory of Theoretical Physics, Joint Institute for Nuclear Research, Dubna 

We examine the large-ir QCD evolution of the twist-three gluonic-pole strength defining 
an effective T-odd Sivers function, where evolution of the T-even transverse-spin DIS 
structure function 32 is multiplicative. The result corresponds to a colour-factor modified 
spin-averaged twist-two evolution. 

Keywords: polarisation; spin; high-energy. 

PACS Nos.: 13.88. +e 

1. Introduction 

1.1. Motivation 

Single-spin asymmetries (SSA's) have long been an enigma in high-energy hadronic 
physics. Before being first observed experimentally, hadronic SSA's were predicted 
to be small for a variety of reasons. Experimentally, however, SSA's turn out to 
be large in many hadronic processes. It was also long held that such asymmetries 
would eventually vanish at large energy and/or p T . Again, however, the SSA's so 
far observed show no signs of high-energy suppression. 

1.2. What is transverse spin? 

Transverse here indicates a spin vector perpendicular to particle momentum (in 
contrast to parallel as in, e.g., the case of the DIS structure function 171). Note that 
the existence of transverse polarisation is itself independent of particle masses — cf. 
the natural (~ 9 %) LEP beam polarisation. The problem of (small) masses only 
arises when seeking measurable transverse-spin effects, which require spin flip. 

1.3. Single-spin asymmetries 

Generically, SSA's reflect correlations of the form s- (p„k), where s is some particle 
polarisation vector, while p and fe are initial/final particle/jet momenta. A typical 
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example might be: p the beam direction, s the target polarisation (transverse with 
respect to p) and fe the final-state particle direction (necessarily out of the p- 
s plane). Polarisations involved in SSA's must therefore be transverse (although 
there are certain special exceptions). It is more convenient to use an helicity basis 
via the transformation 



The appearance of both |+) and |— ) in the numerator indicates the presence of 
a helicity-flip amplitude. The precise form of the numerator indicates interference 
between two helicity amplitudes: one helicity-flip and one non-flip, with a relative 
phase difference, the imaginary phase implying naive T-odd processes. 

Kane et al. [1] realised that in the massless (or high-energy) limit and the Born 
approximation a gauge theory such as QCD cannot furnish either requirement: 
massless fermion helicity is conserved and tree diagram amplitudes are always real. 
This led to the famous statement [1]: ". . . observation of significant polarizations in 
the above reactions would contradict either QCD or its applicability." 

As we now know, large asymmetries were found, but QCD nevertheless survived! 
Efremov and Teryaev [2] discovered one way out within the context of perturbative 
QCD. Consideration of the three-parton correlators involved in, e.g., gi, leads to 
the following crucial observations: the relevant mass scale is not that of the current 
quark, but of the hadron; the pseudo-two-loop nature of the diagrams leads to an 
imaginary part in certain regions of partonic phase space [3] . 

However, it took time before progress was made and the richness of the available 
structure was fully exploited — see [4]. Indeed, it turns out that there are a variety 
of mechanisms to produce SSA's: 

• Transversity: hadron helicity flip is correlated to quark flip. Chirality conserva- 
tion requires another T-odd (distribution or fragmentation) function. 

• Internal quark motion: correlation between the transverse polarisation of a quark 
and its own transverse momentum. This corresponds to the Sivers function [5] 
and requires orbital angular momentum together with soft-gluon exchange. 

• Twist-3 transverse-spin dependent three-parton correlators (cf. gi): the pseudo 
two-loop nature provides effective spin flip and an imaginary part via pole terms. 

The second and third mechanisms turn out to be related. 

2. Single-Spin Asymmetries 
2.1. Single-hadron production 

As a consequence of the multiplicity of underlying mechanisms, there are various 
types of distribution and fragmentation functions that can be active in generating 




(1) 



(2) 



(T|T) + (I|I) (+1+) + (-!-) 



SSA's: 
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• higher-twist distribution and fragmentation functions, 

• fey-dependent distribution and fragmentation functions, 

• interference fragmentation functions, 

• higher-spin functions, e.g., vector-meson fragmentation functions. 

Consider the case of single-hadron production with one initial-state transversely 
polarised hadron: 

A\P A ) + B(P B )^h(P h ) + X, (3) 

where hadron A is transversely polarised while B is not. The unpolarised (or spin- 
less) hadron h is produced at large transverse momentum PhT and PQCD is thus 
applicable. Typically, A and B are protons while h may be a pion or kaon etc. The 
following SSA may be measured: 

A h _ da(S T ) - da (S T ) ,. } 
T da(S T )+da(S T ) [ ' 



AHPa) 



B(Pb) 




Figure 1. Factorisation in single-hadron production with a transversely polarised hadron. 

Assuming standard factorisation to hold, the differential cross-section for such 
a process may be written formally as (cf. Fig. 1) 

d(T = J2Yl Pa'af^a)® fb(x b )®da aa > 77 > ®vfy(z), (5) 

abc aa'^f^f' 

where f a (fb) is the density of parton type a (b) inside hadron A (B), p a aa , is the 
spin density matrix for parton a, Z^/c is the fragmentation matrix for parton c into 
the final hadron h and da aa i^> is the partonic cross-section: 

where M. a p-y8 is the amplitude for the hard partonic process, see Fig. 2. 

For an unpolarised final hadron, the off-diagonal elements of 2?JJ7 C vanish, i.e., 

Tfy^lc ^ ^rr'- Helicity conservation then implies a = a', so there is no dependence 
on the spin of hadron A and all SSA's are zero. To avoid such a conclusion, either 
intrinsic quark transverse motion, or higher-twist effects must be invoked. 
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a, a' 




Figure 2. The hard partonic amplitude, afj-yS are Dirac indices. 



2.2. Intrinsic transverse motion 

Quark intrinsic transverse motion can generate SSA's in three essentially different 
ways (a necessarily T-odd effect): 

(1) for kr in hadron A, f a {x a ) is replaced by V a {x a , kr), which may depend on the 
spin of A (distribution level); 

(2) kt in hadron h allows T> h r j c to be non-diagonal (fragmentation level); 

(3) for k' T in hadron B, fb{xb) is replaced by Vb(xb,k' T ). The spin of b in the 
unpolarised B may couple to the spin of a (distribution level). 

The three corresponding mechanisms are: (1) the Sivers effect [5]; (2) the Collins 
effect [6]; (3) an effect studied by Boer [7] in Drell-Yan processes. Note that all such 
intrinsic-fc^, -«r or ~k' T effects are T-odd; i.e., they require ISI or FSI. Note too 
that when transverse parton motion is included, the QCD factorisation theorem is 
not completely proven, but see [8]. 

Assuming factorisation to be valid, the cross-section is 



The Sivers effect relies on T-odd /cT-dependent distribution functions and predicts 
an SSA of the form 




where again 




E h 



d 3 g(5 , T ) 
d 3 P h 



d 3 a{-S T ) 
d 3 P h 



1^1 ^ / dxadx b Aq f a (x a ,kj.) fb(Xb) 

' * / TV 7. 




da-(x a ,x b ,k T ) 
di 



D h/c (z), (9) 



where Ao 



/ (related to /f^) is a T-odd distribution. 
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2.3. Higher twist 

Efremov and Teryaev [2] showed that non- vanishing SSA's can also be obtained 
in QCD by resorting to higher twist and the so-called gluonic poles in diagrams 
involving qqg correlators. Such asymmetries were later evaluated in the context of 
factorisation by Qiu and Sterman, who studied direct photon production [4] and 
hadron production [9]. This program has been extended by Kanazawa and Koike 
[10] to the chirally-odd contributions. The various possibilities are: 

) <8> fb(xb) ® da ® D h / C (z) 

abc 

+ A T f a (x a ) <g> E b F (x b ,y b ) <g> da ® D h/c (z) 

+ A T f a {x a ) <g> f b (x b ) ® da" ® D$ c (z)} (10) 

The first term is the chirally-even three-parton correlator pole mechanism, as pro- 
posed in [2] and studied in [4, 9]; the second contains transversity and is the chirally- 
odd contribution analysed in [10]; and the third also contains transversity but re- 

(3) 

quires a twist-3 fragmentation function D) , ' , 



2.4. Phenomenology 

Anselmino et al. [11] have compared data with models inspired by the previous pos- 
sible (fcT-dependent) mechanisms and find good descriptions although they cannot 
differentiate between mechanisms. The calculations by Qiu and Sterman [4] (based 
on three-parton correlators) also compare well but are complex. However, the twist-3 
correlators (as in 32) obey constraining relations with fcT-dependent densities and 
also exhibit a novel factorisation property. 



2.5. Pole factorisation 

Efremov and Teryaev [2] noticed that the twist-3 diagrams involving three-parton 
correlators can supply the necessary imaginary part via a pole 
term; spin-flip is implicit (and due to the gluon). The standard 
propagator prescription (— •— in Fig. 3), 




1 



= P-^Ti^(fc 2 ) 



(11) 



Figure 3. Example of 
a propagator pole in a 
three-parton diagram. 



fc 2 ± ie fc 2 

leads to an imaginary contribution for fc 2 — > 0. A gluon with 
x g p inserted into an (initial or final) external line p' sets fc = 
p' — x g p and thus x g — > <^ fc 2 — > 0, see Fig. 4. 

This can be performed systematically for all poles (gluon 
and fermion): on all external legs with all insertions [12]. The 
structures are still complex: for a given correlator there are 
many insertions, leading to different signs and momentum 
dependence. 
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Pole 
Part 



-ITT 

p'.p 



Figure 4. An example of Pole factorisation: p is the incoming proton momentum, p' the outgoing 
hadron and £ is the gluon polarisation vector (lying in the transverse plane). 



The colour structure of the various diagrams (with the types of different 
soft insertions) is also different (we shall examine this question shortly). In all 
cases (examined) it turns out that just one diagram dominates in the large-iV c 
limit, see Fig. 5. All other insertions into external (on-shell) 
legs are relatively suppressed by 1/ A^ 2 . This has been exam- 
ined in detail by Ramilli (Insubria Univ. Masters thesis): the 
leading diagrams provide a good approximation. It has yet to 
be repeated for all the other twist-3 contributions (e.g., also 
in fragmentation). 

A question now arises: what is the relationship between 
the twist-3 and fc^-dependent mechanisms? It might be hoped 
that, via the equations of motion etc., unique predictions 
for single-spin azimuthal asymmetries could be obtained by 
linking the (e.g., Sivers- or Collins-like) /cr-dependent mech- 
anisms to the (Efremov-Teryaev) higher-twist three-parton 
mechanisms. Ma et al. [14] made a first attempt for the Drell- 
Yan process, but the predictions were found not to be unique. Ji et al. [15] have also 
examined the relationships between ^-dependent and higher-twist mechanisms by 
matching in the intermediate kr region of common validity. 




Figure 5. Example of 
a dominant propagator 
pole diagram. 



3. More on Multiparton Correlators 
3.1. Colour modification 

In [16] we provided an a posteriori proof of the relation between twist-3 and the 
Sivers function. The starting point is a factorised formula for the Sivers function: 



dAcr ~ J d 2 k T dx f s (x, k T ) e psPk r Tr[ 7p H(xP, k 7 



(12) 



Expanding the subprocess coefficient function H in powers of k T and keeping the 
first non- vanishing term leads to 



J d 2 k T dx f s (x, k T ) fc£ e psPk ? Tr 



1p 



dH (xP, k T ) 



dk% 



(13) 



fc„=0 
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Exploiting various identities and the fact that there are other momenta involved, 
this can be rearranged into the following form: 



dAa <~ M / dx / W (x) e asPn Tr 



f 



where 



(x) = jc 



dH(xP, k T ) 



2M 2 ' 



fc~=0 



(14) 



(15) 



This coincides with the master formula of [17] for twist-3 gluonic poles in high-p T 
processes. The Sivers function can thus be identified with the gluonic-pole strength 
T(x, x) multiplied by a process-dependent colour factor. 

The sign of the Sivers function depends on which of ISI or FSI is relevant: 

1 



(16) 



where d is a relative colour factor defined with respect to an Abelian subprocess. 
Now, to generate high p T , emission of an extra hard gluon is necessary. Then, 
according to the process under consideration, the FSI may occur before or after 
this emission, again leading to different colour factors. In this sense, factorisation 
is broken in SIDIS, although in a simple and accountable manner. In Fig. 6 wc 
consider the application of this relation to high-p T SIDIS. 






Figure 6. Twist-3 SIDIS tt production via quark and gluon fragmentation. 



3.2. Asymptotic behaviour 

The relation between gluonic poles, e.g., the Sivers function, and T-even transverse- 
spin effects, e.g., gi [18-22], remains unclear. There are model-based estimates and 
approximate sum rules. However, the compatibility of general twist-3 evolution with 
dedicated studies of gluonic-pole evolution (Kang et al. [23], Zhou et al. [24] and 



4th. November 2009 11:12 Ratcliffe 



8 P.G. Ratcliffe and O. Teryaev 

at NLO Vogelsang et al. [25]) is still unproven. In the large- a; limit the evolution 
equations for g 2 diagonalise in the double-moment arguments [26]. For the Sivers 
function and gluonic poles, this is the important kinematical region [4]. 

The gluonic-pole strength T(x), corresponds to a specific matrix element [4]. It 
is also the residue of a general qqg vector correlator bv{x\,x 2 ) [27]: 

rp^ 11+12 ^ 

bv{x\,x 2 ) = h regular part, (17) 

X\ — X 2 

which is defined as 

by( Xl ,X 2 ) = -L J^h e ^-^)+^ e ^"( Pl , 51^(0)^^(^)^2)1^,5). 

(18) 

There is another correlator, projected onto an axial Dirac matrix: 

b A {x u X2) = ^f ^i^e^"^^ 2 (pi, S |V;(0)^7 5 s^(Ai)^(A 2 )b 1 , S ) . 

M J TT 

(19) 

This is required to complete the description of transverse-spin effects, in both SSA's 
and (?2- The two correlators have opposite symmetry properties for x\ <-> x 2 : 

bA(xi,x 2 ) = bA(x2,xi), b v (xi,x 2 ) = -bv(x 2 ,xi), (20) 

which are determined by T invariance. 

In both DIS and in SSA's only a particular combination appears [19]: 

b-(xi,x 2 ) = b A (x 2: xi) - b v (xi,x 2 ). (21) 

The QCD evolution equations [20-22] are usually written in terms of another quan- 
tity, which is expressed as matrix elements of the gluon field strength: 

Y(xi,x 2 ) = (xi - x 2 )b-(xi,x 2 ). (22) 

It should be safe to assume that b-(xi 1 x 2 ) has no double pole and thus 

T(x)=Y(x,x). (23) 

Evolution is easiest studied in Mellin-moment form; for Y(x,y) these become 
double moments: 

Y mn = J dxdyx m y n Y{x,y) 1 (24) 

where the variables are restricted to |y| and \x — y\ < 1. We wish to examine the 
behaviour for x and y both close to unity and therefore close to each other. Thus, 
the gluonic pole provides the dominant contribution: 

YimY{x,y)=T{^) + 0{x-y). (25) 

x, y^l 



In this approximation (now large to, n) the LO evolution equations simplify: 



d r™" = 4(C F + % ) In nY nn , (26) 
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where the evolution variable is s = (3 1 lnlnQ 2 . In terms of T(x) this is 



^>-«(* + %)/>{^(rbr + ™ 



(27) 



which is similar to the unpolarised case, but differs by a colour factor (Cf + Ca/2) 
and a softening factor (1 — z)/(l — x). 

The extra piece in the colour factor (Ca/2) with respect to the unpolarised case 
(Cf) reflects the presence of a third active parton — the gluon. That is, the pole 
structure of three-parton kernels is identical, but the effective colour charge of the 
extra gluon is Ca/2. The softening factor is inessential to the asymptotic solution, 
it simply implies standard evolution for the function f(x) = (1 — x)T(x). For an 
initial distribution /(x,Qq) = (1 — x) a , the asymptotic solution [28] is the same but 
modified by a —> a(s), with 



For T(x), a shifts to a — 1; the evolution modification is identical; the spin-averaged 
asymptotic solutions are thus also valid for T(x). This large-x limit of the evolution 
agrees, bar the colour factor itself, with studies of gluonic-pole evolution [23-25]. 

4. Summary and Conclusions 

Viewing the Sivers function as an effective twist-3 gluonic-pole contribution [16], 
it is seen to be process dependent: besides a sign (ISI vs. FSI), there is a process- 
dependent colour factor. This factor is determined by the colour charge of the initial 
and final partons. It generates the sign difference between SIDIS and Drell-Yan at 
low p T , but in hadronic reactions at high p T it is more complicated. Such a picture 
is complementary to the matching in the region of common validity. Such matching 
between various p T regions now takes the form of a p T -dependent colour factor. It 
also lends some justification to the possibility of global Sivers function fits [29]. 

We have also shown the applicability of generic twist-3 evolution equations to 
the Sivers function. Its effective nature allows us to relate the evolution of T-odd 
(Sivers function) and T-even (gluonic pole) quantities. An important ingredient of 
our approach is the large-x approximation, where gluonic-poles dominate and the 
evolution simplifies. We have found that the Sivers function evolution is multiplica- 
tive and described by a colour-factor modified twist-2 spin-averaged kernel [30]. 
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